Monte Carlo study of gating and selection in potassium channels 
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The study of selection and gating in potassium channels is a very important issue in modern 
biology. Indeed such structures are known in essentially all types of cells in all organisms where 
they play many important functional roles. The mechanism of gating and selection of ionic species is 
not clearly understood. In this paper we study a model in which gating is obtained via an affinity- 
switching selectivity filter. We discuss the dependence of selectivity and efficiency on the cytosolic 
ionic concentration and on the typical pore open state duration. We demonstrate that a simple 
modification of the way in which the selectivity filter is modeled yields larger channel efficiency. 

PACS numbers: 87.10.-c; 87.16.-b; 66.10.-x 



I. INTRODUCTION 

Potassium currents across nerve membranes have been 
longly studied (see, for instance, and the reviews 0- 
8|]). It is now known that ionic channels selecting potas- 
sium currents are present in almost all types of cells in 
all organisms and that they play many important and 
different functional roles. 

An important improvement in experiments was the a 
pearing of the patch clamp technique, see for instance 
which permitted the measurement of the ionic current 
flowing through a single channel on the cell membrane. 
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Different types of measurements, see for instance 
provide a very detailed description of the behavior of 
potassium channels. Less is known on their structure jfj; 
something has been inferred starting from functional be- 
havior and some has been deduced via direct inspection. 

Despite the large variety of ionic channel types, they 
all form selective pores in the cell membrane which open 
and close stochastically and, when in the open state, al- 
low permeation of a selected ionic species (potassium in 
K + -channels) . Gating, i.e., their ability to open and 
close, and selectivity, i.e., their ability to allow the flux 
of a particular ionic species in the cytoplasm, are not yet 
completely understood. The way in which gating and 
selection are achieved can be different from channel to 
channel Q. 

On theoretical grounds this problem has been ap- 
proached with a large variety of methods. Molecular dy- 
namics studies [12| give a very detailed description of the 



phenomenon and their results can also be compared with 



structural experimental information, but they are usually 
not able to provide macroscopic currents estimates due 
to the too small involved time scale. A different approach 
is that of kinetic models 13l4l7| which give very useful 



information since electro-physiological experiments are 
performed over time scales much longer than the atomic 
one. The drawback is represented by the extreme simpli- 
fications on the structure of the channel. 
Models such as the ones in 
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16j describe to some 
extent the dynamics of ion permeation through the se- 
lectivity filter, see also [I?} ■ In this paper, following 18 1, 
we introduce a model where the channel is lumped to a 
two state stochastic point system. A remarkable feature 
of our model is the interaction between the dynamics of 
the ions inside the cell and that of the selectivity filter 
itself. 
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In this spirit, in [18J a very simplified model for the 
pore behavior has been proposed. The main idea is that 
the pore can be in two states, high and Zow-affinity. In 
the former potassium ions can bind to molecules inside 
the pore, while smaller sodium ions cannot. In the latter, 
on the other hand, no ion can bind to the pore. We then 
have that, when the pore is in the low-affinity state, the 
two ionic species permeate in the same way. On the con- 
trary, when the pore is the high-affinity state, sodium 
particles are reflected by the pore, while potassium ions 
bind to it and, eventually, when the pore comes back 
to the low-affinity state, dissociate and possibly exit the 
cell. This selectivity mechanism controls the potassium 
flux through the cell membrane and, hence, yields gating. 
In [rs| it was proven that it is possible to obtain gating 
via selection, but it was also noted that a pronounced re- 
duction in the potassium flux is observed when the frac- 
tion of time the pore spends in the high-affinity state 
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is large; indeed in this case for long period of times no 
permeation is allowed. 

Think for instance of the Kcv from chlorella virus 



PBCV-1 10(, which have no long cytosolic regions able 
to close and open in response to membrane voltage. For 
this kind of channels it is attractive to think that gating 
is indeed realized via selection 18[; in other words re- 



searchers were led to suppose that different current levels 
are indeed obtained through the selection mechanism. 

The basic motivation of this paper is studying the way 
in which the ion dynamics in cytosol affects the pore be- 
havior and exploring the possibility to reduce the potas- 
sium flux loss due to selection. In our model we assume 
that the ionic species diffuse in the cytosol through inde- 
pendent random walks on a two-dimensional lattice and, 
for the potassium and sodium flux ratio, we find results 
comparable to those in |. where the motion of ions 
in cytosol was assumed to be uniform (constant veloc- 
ity) . As already noted in 18[ , good selection is obtained 
when the time fraction spent in the high-affinity state 
is chosen long enough. In this case, the potassium flux 
through the channel is remarkably reduced. Our results 
suggest that, by modelling the cytosolic dynamics as a 
diffusive rather than uniform motion, the potassium flux 
loss is less important. Moreover, we prove that it is pos- 
sible to compensate this effect just allowing the pore to 
stuff more than one particle at a time. For realistic con- 
centrations of potassium inside the cell it is seen that, 
allowing three (in rare occurencese four) potassium ions 
to be packed in the pore, when it is in the high-affinity 
state, is enough to compensate the potassium flux re- 
duction due to the selectivity filter. This hypothesis is 
compatible with what is known about the structure of 
potassium channels 

We also perform an extensive study of the depen- 
dence of fluxes on the two main physiological parame- 
ters, namely, the fraction of time the pore spends in the 
high-affinity state and the concentration of the ions in 
the cell. As it is easy to guess the sodium flux is di- 
rectly proportional to the ionic density in the cell. For 
potassium ions, on the other hand, this linear behavior is 
found only for small values of the time fraction spent in 
the high-affinity state (in this case the two ionic species 
behave similarly), while for large values of such a time 
fraction a sub-linear behavior is found. We studied the 
behavior of fluxes as a function of the time fraction spent 
in the low-affinity state at fixed ionic concentration; we 
found a sub-linear increasing behavior both for sodium 



and potassium. 

Since the pore is modeled in a very simple way (a two 
state random variable independent by the random walk 
dynamics inside the cell), it is possible to study analyt- 
ically the model. We perform an approximate study in 
the two-dimensional case and we are able to explain the 
results found via the Monte Carlo simulation. The most 
interesting effect that is pointed out by this analysis is 
that the sub-linear behavior discussed above is due to a 
sort of depletion of the region close to the pore, which is 
observed when the typical time spent by the pore in the 
low-affinity state is large. The one dimensional model 
has been studied to some extent exactly in order to sup- 
port our approximate two-dimensional analysis and to 
test the reliability of our Monte Carlo computation. 

The paper is organized as follows. In Section |n] we in- 
troduce the model, describe the quantities that are mea- 
sured with the Monte Carlo procedure, and give their 
physical interpretation. Sections IIIIlTVl are devoted to 
the discussion of our numerical results: we first focus 
on the possibility of realizing gating via selection and on 
the way in which the potassium flux loss can be reduced; 
then we discuss the dependence of the fluxes on the phys- 
iological parameters of the model. In Section |VI] we give 
an approximate analytical study of the two-dimensional 
model and in Section [VIII we summarize our conclusions. 
Appendices are devoted to the detailed algorithmic def- 
inition of the model and to the exact one-dimensional 
study. 



II. THE MODEL 



oer is inspired to 



1SJ . We give now 



The model that we study in this par. 
the two-dimensional one proposed in 
a quick glance to the model and postpone the precise 
description of the algorithm to the appendix [A] The 
cytosolic region of the cell is modeled via a finite two- 
dimensional square lattice A with L 2 sites, where L is an 
odd integer (sec figure[T](left)). We consider a two dimen- 
sional model for simplicity, but as it will be discussed in 
the sequel, the model is able to capture the main features 
of the real behavior. Two ionic species perform indepen- 
dent symmetric random walks on the lattice with reflec- 
tivity conditions on the boundary <9A of the volume A. 
One of the site of the boundary is special and it is called 
pore; its behavior, that will be described below, makes 
the potassium walkers not independent. The number of 
potassium and sodium walkers are respectively denoted 
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FIG. 1. On the left: thin line intersections represent the lattice sites; intersections between thin and thick lines represent the 
boundary sites. The pore is in the middle of the right-hand boundary side, (a) Rules for a particle in the bulk: the particle 
jumps with uniform probability 1/4 to one the four nearest neighboring sites, (b) Rules for a particle close to the boundary: 
the particle jumps with uniform probability 1/4 to one the three nearest neighboring sites in the lattice; with probability 1/4 
the particle does not leave the site, (c) Rules for a particle in the corner: the particle jumps with uniform probability 1/4 to 
one the two nearest neighboring sites in the lattice; with probability 1/2 the particle does not leave the site, (d) The rules for 
a particle in the site neighboring the pore depend on the state of the pore and are explained in the picture on the right. Pore 
in the high-affinity state (potassium): when the pore is free the particle jumps with uniform probability 1/4 to one of the four 
nearest neighboring sites (pore included); when the pore is occupied the particle jumps with uniform probability 1/4 to one of 
the three nearest neighboring sites in the lattice (it cannot enter the pore) and with probability 1/4 the particle does not leave 
the site. Pore in the high-affinity state (sodium): the particle jumps with uniform probability 1/4 to one of the three nearest 
neighboring sites in the lattice (it cannot enter the pore) and with probability 1/4 the particle does not leave the site. Pore 
in the low-affinity state: same rule as that for the potassium particle faced to the pore in the high-affinity state; it is worth 
remarking that the pore in the low-affinity state will never be occupied (see the description of the algorithm in the text). 



by Nk and Njq a . The fact that the walkers arc indepen- 
dent on A means that the position of a particle does not 
affect the motion of the others, in particular no constraint 
to the number of particles on each site is prescribed. 

The pore is modeled by the site in the middle of one 
of the four sides of the boundary dA. Two states are 
allowed for the pore: high- affinity and low-affinity. The 
pore switches between the two states randomly; the prob- 
ability that the pore is in the low-affinity state is denoted 
by p. The behavior of the particles on the site neighbor- 
ing the pore depends on the state of the pore itself as it 
is precisely stated in figure [1] The idea is the following: 
if the pore site is in the low-affinity state, both potas- 
sium and sodium ions can enter it; when they enter the 
pore they immediately dissociate so that, with probabil- 



ity 1/2, they reenter A and with the same probability 
exit the system. If the pore is in the high-affinity state, 
sodium ions are reflected by the pore site and potassium 
ions are reflected by the pore if it is occupied, while they 
can enter it if it is free. In this last case the ion does 
not dissociate immediately, but it remains bounded to 
the pore until its state is changed to the low-affinity one. 
As noted above, due to the pore rule, the sodium walkers 
are independent while the potassium ones are not. 

Whenever an ion exits the system a particle of the same 
species is put at random with uniform probability 1/L 2 
on one of the L 2 sites in A so that the number of sodium 
and potassium ions in the system remains constant. 

The model described above is implemented with the 
Monte Carlo algorithm which will be described in detail 
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in the appendix [3] An iteration or sweep of the Monte 
Carlo is the collection of the steps that are performed at 
each time t. The number of iterations that are performed 
in a simulation is t max . 



A. Measured quantities and goal 

The particles that exit the system model the outgoing 
ionic flux. This is precisely what we want to measure. 
For a single numerical experiment we let M^ a (t) (resp. 
Mk(£)) be the number of sodium (resp. potassium) par- 
ticles which exited the system in the time interval [0, i]. 
Moreover, we let the sodium (resp. potassium) flux /Na(t) 
(resp. /k(£)) at time t be M^ a (t)/t (resp. M^{t)/t). Since 
we defined the Monte Carlo algorithm in such a way that 
the number of sodium and potassium ions keep constant 
in the volume A, both f^ a (t) and /k(£) approach a con- 
stant value for t — > oo. To estimate these limiting val- 
ues, which are denoted by /N a and /k, we perform dif- 
ferent numerical experiments with number of iterations 
tmax sufficiently large, we then average the experimen- 
tal measures /Na(imax) and /K(imax) and compute the 
associated statistical errors. 

Our aim is to measure /Na and /k and to prove that 
the model provides selection in an efficient way. The two 
ionic species interact in the same way with the pore in 
the low-affinity state, see the rules in figure [T] hence, 
when p is close to one, no difference should be observed 
in the potassium and sodium fluxes. A different behavior 
should be observed at smaller p, indeed when the pore is 
in the high-affinity state the sodium flux is blocked while 
some residual potassium flux should be recorded. 

We shall then compute the ratio ik/ika and prove 
that it becomes large for p small. Moreover we shall 
discuss if this selection mechanism is efficient, namely, 
we shall see that selection is provided with a potassium 
flux comparable with the flux measured when the pore is 
constantly in the low-affinity state. In particular in the 
sequel we shall discuss a slightly different model which 
will provide selection without flux reduction. 



B. Choice of the parameters 

We first describe how the physiological parameters of the 
model have been chosen. We have considered different 
values of the number of potassium (Ak) and sodium 
(A^Na) particles. Although cytosol is typically much 



richer in potassium than in sodium, in order to check the 
efficiency of the selectivity filter, we have always taken 
ANa = Ak- In all the simulations discussed in this sec- 
tion we have chosen L = 10 1, while we considered the 
cases 



A Na = A K = 100, 1000, 3000, 5000, 10000 



and 



p = 0.001, 0.003, 0.005, 0.007, 0.01, 0.03, 
0.05,0.07,0.1,0.3,0.5,0.7,1 

Note that the total time fraction the pore spends in the 
low-affinity state is essentially equal to t max p. 

It is worth noting that, provided L = 101, the most 
physiologically reasonable choice for the number of potas- 
sium ions among those listed above seems to be Ak = 
100. This is supported by the following very rough com- 
putation. The concentration of potassium ions in cytosol 
is approximatively equal to c = 150 mM (millimolar) , 
that is c mole per cubic meter. The radius of the potas- 
sium ion is r = 1.52 x 10 -10 m; assuming that each site 
of the lattice can accommodate at most n ions, we can 
ascribe the volume n(2r) 3 to each site of the lattice. The 
number of potassium ions that must be put on the lattice 
is then 

number of potassium ions — cL 2 n(2r) 3 n,A = 25.88n 

where tia = 6.02 x 10 23 is the Avogadro number, i.e., 
the number of ions in a mole. Assuming that reasonable 
values of n ranges between 1 and 10, we have that realistic 
values of the number of potassium ions on the lattices 
ranges from 25 to 250. 

The sole Monte Carlo parameter to be fixed is the num- 
ber of sweeps t max . As explained above it must be taken 
large enough in order to get a good estimate of the out- 
going fluxes. In all the cases that we considered it was 
sufficient to take i max = 10 7 . 



III. GATING AND SELECTION 

As it has already been stated our first purpose is check- 
ing the behavior of the ratio /k/ /Na as a function of the 
probability p. Results are shown in figure [2j from the 
bottom to the top the ratio /k/ /Na for AN a = Ak 
11)000. 5000, 3000, IOOO, IOO has been plotted as a func- 
tion of p. We have that the ratio between the potassium 
and the sodium outgoing fluxes increases when the prob- 
ability for the pore to be in the low-affinity state is de- 
creased. This is precisely the result in 18(; in particular 
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FIG. 2. Log-log plot of the ratio /k//ni as a function of 
the probability p. Data refers to the cases AWa = Nk = 
10000(o), 5000(O), 3000(A), 1000(x), 100(+). 

for A^Na = Nk = 100 our data are not only qualitatively 
but also quantitatively close to those in Jl8| , although 
a different ion dynamics in the cytoplasm has been con- 
sidered. As noted in [l^ it is then possible to imagine 
a potassium channel in which gating is realized via the 
selectivity filter. 

The log-log plots in figure [2] shows a straight line be- 
havior for p large, say p > 0.08. Note that the threshold 
where the linear behavior starts is smaller for smaller val- 
ues of the ionic density in cytosol. This means that in 
this regime, say p > 0.08, the ratio /k//n3 can be well 
approximated by a power law 

i^ = a 1 {pr bi a) 

jNa 

The coefficient ax and b\ can be computed by 
fitting the data plotted in figure [2J we obtain 
a x = 0.863,0.904,0.974,1.014,1.072 and bx = 
0.860,0.798,0.694,0.616,0.480 respectively for N K = 
A^Na = 100, 1000, 3000, 5000, 10000. 

It is worth noting that, if the density of ions in the 
cytosolic region is not too high, good selection ratios are 
reached at not too small p; this is good news, since too 
low values of p would imply very small ionic currents. In 
particular, good results are found for the physiologically 
reasonable parameters A^Na = Nk = 100. When the 
number of particles in the system is raised the selection 
effects are widely reduced. This is quite intuitive, indeed, 
when the density is large, the average number of particles 
trying to enter the pore is large, as well; in other words 
the time between two consecutive attempts to enter the 



FIG. 3. Log-log plot of the normalized fluxes as func- 
tions of the low-affinity state probability p. For the 
potassium the data in the picture refers to the cases 
N K = 10000(o), 5000(O), 3000(A), 1000(x), 100(+); for the 
sodium we used the symbol □ for the five cases A<Na = 
10000, 5000, 3000, 1000, 100. 

pore is small. If the typical duration of the intervals in 
which the pore stays in the high-affinity state is larger 
than this "knocking" time, a lot of attempts to enter the 
pore will abort since the filter is occupied by another 
particle. It is then clear that the selection mechanisms 
does not work in the optimal regime. In the sequel we 
shall propose a modified version of the model which will 
take care of this problem. 

We have seen that the algorithm provides good selec- 
tion for p small enough; at physiological concentrations, 
in order to get a ratio /x/ZNa larger than 100, the low- 
affinity state probability p must be chosen smaller than 
0.01. In this regime we expect a potassium flux reduction 
with respect to the case p = 1. We estimate, now, this 
loss by plotting in figure [3] the normalized potassium and 
sodium fluxes, namely, the ratio between fluxes at any p 
with the corresponding fluxes at p = 1. We first note that 
the flux reduction is more effective for the sodium ions 
rather than for the potassium ones; this is obvious, since 
the selectivity filter favors potassium particles. More- 
over, in the case of sodium particles, the reduction does 
not depend on the number of particles on the lattice and 
this is due to the fact that all the sodium related mea- 
sured quantities are directly proportional to the sodium 
density in cytosol. 

Potassium flux, on the other hand, experiences a re- 
markable loss at small low-affinity state probability p. 
This reduction gets more and more important as the 
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density is increased. At physiological potassium den- 
sity the flux reduction is about the 40% at very small 
low-affinity state probability. In respect to this point it 
seems that our algorithm is more efficient than the one 
proposed in 18[, indeed, in that paper, with a similar 
selection ability, larger potassium flux losses were found. 
In other words the importance of potassium flux loss de- 
pends on how the cytosol dynamics is modeled, indeed in 
our model ions move inside the cell according to a sym- 
metric random walk, while in [hsl ] they perform a uniform 
motion until the boundary is reached. In the next sec- 
tion we propose a modified version of the model aiming 
to reduce potassium losses. 
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IV. IMPROVING EFFICIENCY 

As it has been discussed above, the model which has been 
proposed before is able to describe the gating attitude of 
a potassium channel via a suitable selectivity filter. 

The main problem with this approach is that the pres- 
ence of the filter reduces the amount of the permeating 
potassium ions. To be more precise, in order to obtain 
a remarkable selection one has to assume that the time 
fraction in which the pore is in the low-affinity state is 
very small. From the data in figure [2l it follows that, 
at physiological concentrations, in order to get the ratio 
,/k//Na larger or equal to 100 one has to assume p < 0.1. 
As shown in figure El at these values of the low-affinity 
state probability an important reduction of the potas- 
sium flux is observed. 

This phenomenon has a simple explanation: at small 
p many attempts performed by K + ions to exit the cell 
abort since it often happens that the pore is in the high- 
affinity state and occupied by another potassium particle. 
This problem can be bypassed by assuming that more 
than one ion at time can be accommodated in the pore. 
We test this idea in the simplest possible way: we do 
not model the pore structure or any possible interaction 
between ions in the pore as it has been done in 18[ (this 
is obviously a very interesting question deserving future 
investigation). We simply assume that there is no upper 
bound to the number of ions that can be accommodated 
in the pore (see below). In other words when the filter is 
in the high-affinity state, a potassium particle can enter 
it even if the pore is currently occupied by one or more 
potassium ions. 

This modified version of our model will assure no 
potassium flux reduction since all the potassium ions try- 



FIG. 4. Log-log plot of the ratio /k//ns as a function of the 
probability p for the modified model; this picture should be 
compared with figure[2] Data refers to the cases A^Na = Nk = 
10000(o), 5000(O), 3000(A), 1000(x), 100(+). 



ing to enter the pore in the high-affinity state will be 
accepted and then will be released when the state of the 
filter will switch to the low affinity one. 

This model will not be physiologically reasonable if 
the typical number of potassium ions stuffed inside the 
pore is large. From what it is know n on th e structure of 
ionic channels (see, for instance 
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19|, 2(3]), a potassium 



channel is about 12 A long and 2.5 A in diameter. This 
means, recalling that the radius of a potassium ion is ap- 
proximatively 1.5 A, that at most four or five potassium 
particles can coexist inside the pore. We stress that in 
this simple model we do not take into account any inter- 
action, for instance electrostatic repulsion, between ions 
in the pore. As we shall see, at physiological densities, 
i.e., ANa = Ak = 100 and for p > 0.01, the average 
number of potassium ions in the pore will be very low 
and at most four particles (in very rare cases) will be 
simultaneously accommodated inside the pore. We shall 
then conclude that in this regime our simple model is 
reasonable on physiological grounds. 

The results are shown in figures |4] and [5] By comparing 
figures |3] and [21 it is clear that the selectivity attitude of 
the new model does not depend on the density of the ions 
in cytosol and that the modified model is more selective. 
Figure [5] is even more interesting, indeed it shows that 
the potassium flux does not depend on the low-affinity 
state probability. In other words the model introduced in 
this section behaves precisely as we expected and, hence, 
it provides the solution to the potassium flux loss. 
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FIG. 5. Log-log plot of the normalized fluxes as func- 
tions of the low-affinity state probability p for the modi- 
fied model; this picture should be compared with figure [3] 
For the potassium the data in the picture refers to the cases 
Nk = 10000(o), 5000(O), 3000(A), 1000(x), f00(+); for the 
sodium we used the symbol □ for the five cases AWa = 
10000,5000,3000,1000,100. 



It is important, now, to discuss the physiological rea- 
sonableness of the model. In order to answer this ques- 
tion, which is very important in the present context, 
we have run a long simulation (about 3 x 10 7 sweeps) 
recording at each instant the number of potassium par- 
ticles accommodated in the pore site. In the three cases 
A^a = N K = 100,1000,5000 (see figure El (a), (b), and 
(c)) the number of simultaneously stuffed particles is, 
along the whole simulation, respectively smaller than 4, 
15, and 62. 

It then follows that at physiological values of the cy- 
tosolic ion density, i.e., for N-^ a — — 100 and for 
low-affinity state probability not too small (p > 0.01), 
the model is completely reasonable. 

We can then conclude that a model of potassium chan- 
nels is feasible in which the gating attitude is obtained 
via a selectivity filter and such that the filter does not 
produce any potassium flux reduction. 



V. DISCUSSION 

We now discuss the behavior of the sodium and potas- 
sium flux in the basic model as function of the two main 
physical parameters: the ionic cytosolic density and the 
time fraction in which the pore is in the low-affinity state 
(which, we recall, is essentially p). In this section we 
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TABLE I. Values of the parameters 02 and 62 introduced in 
equation ([2]) and obtained by fitting the data plotted in fig- 
ure [3 

shall simply describe the numerical results; the next sec- 
tion will be devoted to their physical interpretation via 
a simple analytical model. 

It is difficult to have an intuition on the potassium flux 
behavior, indeed the way in which the selectivity filter 
acts upon its current is not transparent at all. Since 
the filter acts simply as a gate on sodium, it is rather 
natural to expect that the sodium flux will be directly 
proportional to the sodium density on the lattice; on the 
other hand the way in which it depends on the the low- 
affinity state probability, i.e., on the time fraction the 
filter spends in the low-affinity state, cannot be easily 
guessed a priori. 

We discuss first how the sodium flux depends on con- 
centration in cytosol. In figure [7] we have plotted the 
sodium flux /n a as a function of the sodium density 
0Na = N^ a /L 2 on the lattice. Only the curves corre- 
sponding to the values p = 1,0.7,0.5,0.3,0.1 are shown; 
other values of the low-affinity state probability yield 
similar results. As expected, the flux is directly propor- 
tional to the density with slope depending on the low- 
affinity state probability. We have fitted our data with 
the function 



/Na — 0-2 + 62gNa 



(2) 



and computed the parameters 02 and 62 for the different 
values of p; results are shown in the table HI It is very 
interesting to remark that the slope coefficient 62 depends 
linearly on p in the interval [0, 0.1]; for larger values of the 
low-affinity state probability the behavior departs from 
linearity. To illustrate this remark, in figure [8] we have 
plotted the coefficient 62 as a function of p (we recall that 
the time fraction spent by the pore in the low-affinity 
state is, roughly speaking, equal to p). The picture shows 
that for small low-affinity state probabilities the slope 62 
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(a) (b) (c) 

FIG. 6. Number of potassium ions accommodated inside the pore in an interval of time of few thousands sweeps out of3xl0 7 for 
the modified model. Parameters of the simulation: p = 0.01, iV N a = N K = 100 (a), iV N a = N K = 1000 (b), iV N a = N K = 5000 
(c). 
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FIG. 7. Sodium flux /n 3 as a function of the sodium den- 
sity £Na- Data in the picture refers to the cases p = 
l(o), 0.7(O), 0.5(A), 0.3(x),0.1(+). The other values of p are 
not shown: the behavior is linear with smaller and smaller 
slope (see table [J). Note that the behavior in the picture is 
much reasonable since the selectivity filter does not affect the 
sodium dynamics. 

is directly proportional to p, that is to the time fraction 
that the pore spends in the low-affinity state. This result 
is consistent with a theorem proven in 21[ in a three- 
dimensional diffusion model. 

The predictions of ([2|) with the coefficient a 2 , 62 chosen 
as in the caption of figure [51 that is 

P 



fx 



(3) 



4 + b 'sP 

are compared with the experimental data in figure |H] and 
the matching is striking. 



FIG. 8. Plot of the slope coefficient 62 introduced in equation 
J5J as a function of the low-affinity state probability p. The 
dotted straight line has been obtained by fitting the data with 
the function 62 = 03 + 63P on the interval [0,0.1]; the result 
of the fit is as — and 63 = 0.118. The solid line has been 
obtained by fitting the data with the function 62 = p/(o-3 + 
b'sp) on the whole data set, i.e., on the interval [0, 1]; the result 
of the fit is a' 3 = 7.96 and 63 = 5.98; note that the derivative 
in is I/03 = 0.1256 w & 3 - 



We now come to the results related to the potassium 
flux. Potassium flux as a function of potassium density 
in the lattice bulk §k = Nk/L 2 has been plotted in fig- 
ure [TDJ In this case the behavior is approximatively lin- 
ear for values of the low-affinity state probability p close 
to one. But at smaller values of p the graph remark- 
ably departs from linearity. In other words, when the 
time fraction spent by the pore in the low-affinity state 
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FIG. 9. Log-log plot of the sodium flux /n s as a function of FIG. 11. Log-log plot of the potassium flux /k as a function 

the low-affinity state probability p. Data refers to the cases of the low-affinity state probability p. Data refers to the cases 

iV N a = 100(»), 1000(x), 3000(A), 5000(D), 10000(+). Lines iV K = 100(»), 1000(x), 3000(A), 5000(D), 10000(+). 
are the graph of the function Q with the appropriate values 

Of £>Na- 

VI. ANALYTICAL STUDY 
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FIG. 10. Potassium flux /k as a function of the potas- 
sium density qk- Data in the picture refers to the cases 
p = 0.005(0), 0.05(»), 0.07(D), 0.1(+), 0.3(x), 1.0(o). The 
other values of p are not shown: the behavior is similar. Lines 
are eye guides. 



is large, the potassium flux is, with good approximation, 
directly proportional to its density in the bulk. On the 
other hand, when the time fraction spent by the filter 
in the low-affinity state is smaller, the potassium flux 
exhibits a sub-linear behavior with density. 

Finally, in figure [TT] the potassium flux as a function 
of the low affinity state probability p is reported. 



In this section we study analytically the basic model in- 
troduced in Section|TI]in order to give a physical interpre- 
tation of the behavior of sodium and potassium flux on 
the species density and on the low-affinity state proba- 
bility. In our model the pore is modeled in a very simple 
fashion, indeed it is just a two state Bernoulli process; 
the main difficulty is, obviously, the interaction between 
the random walk inside the volume A and the pore itself. 
We consider the stationary state reached by a walker and 
denote by gx (resp. gNa) the probability for a potassium 
(resp. sodium) ion to occupy the site x of A neighboring 
the pore. 

Since sodium particles can exit the system only when 
the pore is in the low-affinity state, the sodium flux is 
simply given by 



/Na = o^Na^NaP 



(4) 



where, we recall, A^ a is the number of sodium ions in the 
volume A. Indeed, the probability for a sodium ion to be 
on the site x neighboring the pore is gNa, the probability 
for such a particle to jump to the pore itself is 1/4 and, 
once in the pore, the probability to really exiting the 
system is 1/2. 

When dealing with potassium, one has to take into 
account that, when the pore switches to the low-affinity 
state, the possibly trapped potassium ion is released with 
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FIG. 12. Sodium flux as function of the low-affinity state 
probability. Pluses are the Monte Carlo results for JVns = 
3000. The solid, long dashed, and short dashed lines are 
the corresponding graphs of the functions (p) , (p) and 
/iNa (?) respectively. 

probability one half. We then have for the potassium flux 



1 



1 



where we have denoted by r the stationary probability 
for the pore to be occupied by one of the Nk potassium 
particles. 

The problem, now, is that of computing the station- 
ary probabilities gNa,<?K, and r. The following study is 
related to the two-dimensional model discussed in the 
above sections; in the appendix [B] we address the same 
problem in the one dimensional case, which is an interest- 
ing case, since some of the involved computations can be 
performed exactly. The related results will both support 
our two-dimensional approximate discussion and confirm 
the reliability of our Monte Carlo measurements. 

We first consider the sodium case. An obvious guess is 
that the probability (?Na is = 1/|A|, i.e., it is uniform 
throughout the system. In this way we get 



* (°) 1 at 1 



(6) 



Note that this simple model provides the estimate 8 for 
the constant a' 3 — 7.96, obtained by fitting the data in 
figure [3 and for the constant b' 3 = 5.98. The compar- 
ison between the Monte Carlo data and the theoretical 
prediction is good for small values of the low-affinity 
state probability, while (see figure [T2"]) , the experimental 
data departs from the predicted behavior for large p. 



This behavior is not surprising: at large p the outgo- 
ing flux is so important that one has to expect a sort of 
depletion of the region close to the pore. To estimate 
this effect we develop a sort of "mean field" argument: 
we imagine that in all the sites of A different from the 
one neighboring the pore the stationary probability qj<!a.,o 
is uniform. To get the stationary probability q^ a in the 
site x close to the pore we make the following probability 
balance 



1 



<ZNa = -</Na,0 



(1 -_p)<?Na 



1 1 

42 



P<?Na 



1 1 

42L 2 



PqNa 



where, on the right hand side, the first term is the prob- 
ability that the particle moves from A\{a;} to x, whereas 
the other terms take into account the probability that the 
particle in x remains in x. The above equation, together 
with the obvious normalization condition, 



(L 2 - l)<?Na,0 + 9Na = 1 



yields 



<?Na 



1 

L 2 +p(L 2 -1) 2 /(6L 2 ) 



(5) for the occupation probability and, hence, 



8 i,Na L 2 +p(L 2 -l) 2 /(6L 2 ) 
for the sodium outgoing flux. This approximation pro- 



ides the estimate 8 for the constant 



7.96, obtained 



fxl; this is reasonable, 



by fitting the data in figure U and 4(L 2 - 1) 2 /(3L 4 ) = 
1.33 for the constant b' 3 = 5.98. The graph of the oc- 
cupation probability is plotted in figure 1131 the guessed 
depiction effect is found and the occupation probability 
decreases when p is increased. 

For p small and L large, 
since in this limit the depletion effect is negligible. As 
it is shown in figure [12] the function is a quite good 
approximation of the experimental data, which, as we 
have seen above, are perfectly fitted by the function /N a 
given in equation ([3]). 

A more detailed description of the dynamics close to 
the pore can yield a better estimate of the depletion ef- 
fect. Consider one of the N^ a walkers on A and let, as 
we did before, (?Na be the probability that it occupies the 
site x in A, neighboring the pore. Say that the boundary 
is along the north-south direction and that the pore is to 
the east of the site x. Let <?Na,i be the probability that 
the particle is on the site to the north of x and that it 
is equal to the probability that the walker occupies the 
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site to the south. Let be the probability that the 

particle is on the site to the west of x. Assume that for 
the remaining L 2 — 4 sites of the lattice the occupation 
probability is <7Na,u- 

As before it is not difficult to write the following prob- 
ability balance equations 

<?Na,l = 2g Na ,o/4 + <?Na/4 + <7Na,l/4 
<?Na,2 = 3g N a,o/4 + <?Na/4 

<?Na = 9Na,2/4 + 2(7 Na ,i/4 + (1 - p)<7Na/4 + P<?Na/8 

where we have neglected the terms proportional to 1/L 2 . 
The above equations, together with the obvious normal- 
ization condition 

(L 2 - 4)g Na , + <?Na,2 + 2q N a,l + 9Na = 1 



yield 



(2) 
<?Na 



12/23 



1 + 12(1 + 6p/25)(L 2 - 23/12)/23 



for the occupation probability and, hence, 



- 1 N 



1 



6p(L 2 - 23/12)/25* 



for the sodium outgoing flux. This approximation pro- 



ides the estimate 8 for the constant 



7.96, obtained 



by fitting the data in figure^ and 48(l-23/(12L 2 ))/2E 
1.92 for the constant b' 3 = 5.98. Note that for p small and 
L lar g e > /] 

f(2) 



(2) 
Na 



AO) 

JNa ' 



As it is shown in figure [12] the func- 
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FIG. 13. Referring to the left side of the picture, from the 
top to the bottom the graphs of the functions q^l(p), 9n2(p)i 
q^(p) and (p) for N K = 10000, qg> (p) and q^{p) for 
iV K = 100, are plotted. 
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tion /j^ a is a better approximation of the experimental B 
data with respect to and /^j? . 

The graph of the occupation probability (p) is plot- 
ted in figure [13] the depletion effect is more pronounced 
with respect to that described by the approximation 

A)- 

We now come to the computation of the potassium 
flux. We consider a potassium walker and we note that 
for such a particle the random walk is performed in 
A U {pore}. Wc assume that in the stationary state the 
probability that the potassium walker is in the pore is 
r/Nji- Thus in the stationary state 

^ = jd - P )d-r) qK + 

where, wc recall, qK is the probability that the potassium 
ion occupies the site a; in A neighboring the pore. Indeed, 
by stationarity, the probability that the particle is in the 
pore is equal to the probability that the pore is empty in 
the high-affinity state and that the particle jumps from 
x to it plus the probability that the pore is occupied and 
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FIG. 14. Potassium flux as function of the low-affinity state 
probability. Pluses are the Monte Carlo results for Nk = 
3000. The solid, broken (large), broken (small), dotted lines 
are the corresponding graphs of the functions f^\p), /k^(p)> 
and /k 2 '(p) respectively. 

in the high affinity state. From the equation above we 
easily get 



1 



-ZVk<2k(i -p)- 



(7) 



4~' IV " IVV I ' J p + N K q K (l-p)/4 

Then we have to estimate q^. As we have done in the 
case of sodium we first consider the stationary occupa- 
tion probability uniform inside the volume A and equal 
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to q^p = 1/L 2 . The corresponding probability r 1 - * 1 for 
the pore to be occupied is given by the equation ([7]) and 
the corresponding potassium flux \ given by (0, is 
plotted in figure lMl The comparison with the experimen- 
tal data is very good at small low-affinity probability p, 
while the experimental results depart from the predicted 
behavior at large p. 

As before this is due to the depletion of the region close 
to the pore. We repeat the same computation described 
in detail for the sodium particles. Here we just give the 
main formulas: the first approximation consists in the 
probability balance equation 



<7K = ^<7K,o + g<7KP + ^IKP + 4<?k(1 - P)r 
1 r 1 r 



By using the obvious normalization condition (L 2 — 
1)<7k,o + qK + t/Nk = 1, by neglecting all the terms pro- 
portional to 1/L 2 and to 1/Nk, and by recalling equation 
0, we get 



(1) _ AD + Bp - ^/{AD + Bp) 2 - AACDp 



and 



'Ik 



(i) 



2C P 



A 



B - CrW 



where 



.4 



3 1 



4L 2 



-, B = A + l--p, C = -(l- P ) 



and 



D = CN K 



The corresponding potassium flux fjp given by ([5]) is 
depicted in figure 1141 The improvement with respect to 
is striking. 

We can further improve the way in which the depletion 
effect is estimated as we did for the sodium particles. By 
neglecting, as we did above, contributions proportional 
to 1/L 2 and to 1/Nk, we get the probability balance 
equations 

2gK,o/4 + g K /4-3<? K> i/4 = 
3<2Ko/4 ~ QK,2 + qn/4: = 
g K , 2 /4 + 2<j K ,i/4 - 3g K /4 - q K p/$> 
-(l-p) gK (l_r)/4 = 

and the obvious normalization condition 







(L 2 - 4)q K ,o + Qk. 2 + 2<7k,i + 9k = 1 



The above equations, together with (JT)), yields for 
and q^P the same expression as those for and q^ 
with 

„ 25 1 , 23 „ 37 1 

48 L 2 - 23/12 and B = T2^ + 48"8 P 

whereas C and D are as before. The corresponding potas- 
sium flux , given by ([5]) and plotted in figure 1 141 im- 
proves the estimate 

We finally comment on the behavior of the occupation 
probability as function of the low-affinity state prob- 
ability p. In figure [13] we have plotted q^ and q^ for 
N K = 100 and iV K = 10000. We note that those curves 
tend to the same value of the corresponding curves for the 
sodium occupation probability when p tends to one. This 
proves that the approximations are coherent, indeed, in 
this limit sodium and potassium particles behave simi- 
larly. We also note that sodium occupation probabilities 
are decreasing functions of the low-affinity state proba- 
bility p. For potassium ions, the occupation probability 
decreases at small p and increases at large p. This is due 
to the fact that at large p the probability r that the pore 
is occupied is small and so is the associated contribution 
to the potassium flux. 

VII. CONCLUSIONS 

Potassium channel are special transmembrane proteins 
allowing selected potassium permeation outside cells. 
Gating and selectivity attitudes are not fully understood. 
For special K + -channel types it has been supposed that 
gating is realized via a selectivity filter. In this paper we 
have proposed and studied a lattice model, in the same 
spirit of 1^, in which a pore selection rule ensures both 
selection and gating. This model has been studied both 
by means of a Monte Carlo method and via an analytical 
approximation. We have estimated both potassium and 
sodium fluxes through the membrane. 

We have proven the possibility to achieve selection and 
gating and we have discussed the potassium flux reduc- 
tion due to the presence of the selectivity filter. In partic- 
ular we have shown that, allowing more than one potas- 
sium ion at time to be accommodated inside the pore, the 
potassium flux loss is completely recovered when physi- 
ological cytosolic densities are considered. We then con- 
clude that gating can be achieved via a selectivity filter 
in an efficient way. 

We have also studied extensively the model discussing 
the behavior of potassium and sodium fluxes as functions 
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of the interesting physiological parameters, i.e., the ionic 
cytosolic density and the time fraction the filter spends 
in the low-affinity state. It is remarkable to notice that 
the flux of the ionic species which is not affected by the 
selectivity filter, sodium in our model, is directly propor- 
tional to the cytosolic density. A sub-linear behavior is 
found for potassium 

In conclusion we think that lattice models provide es- 
sential and tunable tools to investigate efficiently the 
properties of ionic channels by analyzing the wealth of 
experimental data available on this subject. 

Appendix A: Detailed definition of the model 

In this section we describe in detail the Monte Carlo 
scheme that we have studied and whose behavior has 
been discussed above. 

We consider an integer time variable t. We set t = 
and choose at random with uniform probability 1/L 2 
the position of the N^ a sodium particles and the Nk 
potassium particles. We then repeat the following steps 
until t equals the given integer number t max : 

1 . set t equal t + 1 ; 

2. select at random the state of the pore: choose the 
low-affinity state with probability p and the high 
affinity one with probability 1 — p; 

(a) if the pore is in the low-affinity state and it is 
occupied by a particle, the particle is released 
with the following rule: it jumps with proba- 
bility 1/2 to the site of A neighboring the pore 
or (with the same probability) it exits the sys- 
tem; 

(b) if a particle exits the system, a particle of the 
same species is put at random with uniform 
probability 1/L 2 on one of the L 2 sites in A; 

3. the position of each particle on the lattice is up- 
dated following the rules defined in figure [TJ 

(a) if a particle enters the pore and the pore is in 
the low-affinity state, the particle is immedi- 
ately released by the pore with the following 
rule: it jumps with probability 1/2 to the site 
of A neighboring the pore or (with the same 
probability) it exits the system; 

(b) if a particle exits the system, a particle of the 
same species is put at random with uniform 
probability 1/L 2 on one of the L 2 sites in A. 



We considered a second model with the aim of im- 
proving the potassium flux efficiency at low-affinity state 
probability. This model is defined as the basic one with 
a single difference: the rule for a potassium ion in the 
site neighboring the pore when the pore is occupied and 
in the high-affinity state (see the right top part of fig- 
ure [T]) is the same experienced when the pore is free in 
the high-affinity state. 

Appendix B: One dimensional model 

The basic model introduced in section Ql] can be eas- 
ily specialized to the one-dimensional case. The lattice 
is A = {1, . . . , L} and the pore is the site L + 1 of Z. 
Equations (j4j and (J5]) become 

ha. = ^A Na g Na p (Bl) 

and 

,/k = \ N K qkP + \pr (B2) 

In order to compute (?Na we consider the random walk 
of a single sodium ion on the lattice {1, . . . , L}. The 
transition matrix elements different from zero are 

ttm+i = \ for i = 1 
for nearest-neighbor jumps to the right, 

7T;+:u = \ for i = 1, . . . , L - 2 
for nearest-neighbor jumps to the left, 

for the nearest-neighbor jump to the left starting from 
the site L, 

n L ,i = -^-j-p for % = 1, . . . , L — 2 

for left jumps starting from the site L (this not zero tran- 
sition matrix elements are due to the fact that particles 
exiting the system are put at random with uniform prob- 
ability in one site of the lattice A), and finally 

1 , 111, 

7Ti,i=2 7rL ' L = 4 P + IZ P + r P ' 

We let qi, with i = 1, . . . , L, be the stationary proba- 
bility that the walker occupies the site i. Since the pore 
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is the site L + 1 of Z, we have that gN a is nothing but g^. 
Recalling the definition of stationary measure 



L 



for i = l, 



, L, we have 



92 = 



2?2 

1 

2 q3 



pqL 



l 

4L 
IL ML 



(B3) 



1 



qL-i 



41 



1 



111 

<7l = 2 qL - 1 + 1 pqL + 4L PqL + 2 ^ ~ 



By using the first L — 1 equations (|B3[) by recursion we 
get 



g fc = <7fc+i 



21 



(B4) 



for k = 1, . . . , L — 1. By recursion, again, it is also easy 
to prove that (|B4j) implies 



QL-i 



2iL - i(i + 1) 



(B5) 



for i = 1, . . . , L — 1, so that the probabilities </i, . . . , qL—i 
are all expressed in terms of q^. Note that qL-i de- 
creases when i runs from L — 1 to 1, so that the depletion 
phenomenon discussed in section IVT1 is found also in the 
one-dimensional case. 

In order to find qx, we finally require that the normal- 
ization condition 



<h 



q L = i 



is satisfied. A straightforward computation yields 

1 



L+(2L 2 - 3£+l)p/12 



(B6) 



The one-dimensional sodium flux can be finally com- 
puted by using (|B1|) with gNa given by (|B6j) . Monte Carlo 
results are compared with this theoretical prediction in 
figure n~5l and the matching is perfect. 

The computation of gx is more difficult, since the inter- 
action between the potassium ions and the pore is not as 
trivial as for sodium. In particular, due to the possibility 
that a potassium ion is trapped in the pore, the potas- 
sium walkers cannot be considered independent. We shall 
treat this case by assuming that in the stationary state 
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FIG. 15. Sodium flux as function of the low-affinity state 
probability in the one-dimensional case for L — 100. The 
points *, x , and + are the Monte Carlo results respectively for 
iV N a = 3000, 1000, 100. The three lines are the graphs of the 
function (|B1|) with gNa given by (|B6[) for the corresponding 
values of iVNa. 



the probability for the pore to be occupied by one of 
the Njt potassium ions is r. We then study a potassium 
walker on the lattice A = {1, . . . , L, L + l}, where the site 
L + 1 is the pore, and at the end we shall require that 
g^+i = r/Nj£- In other words we shall finally assume 
that, in the stationary state, the probability for a single 
walker to occupy the pore is equal to the probability that 
the pore is occupied by one of the walkers divided by the 
number of potassium ions. 

We do not list all the not zero transition matrix ele- 
ments, but only those differing from the sodium case: 

ttl+M = TjjV for i = 1, . . . , L - 1 

due to particles released by the pore when it flips to the 
low-affinity state, and 

= p/4 + p/(4L) + (1 - p)r/2 
kl.l+i = (l-j>)(l-r)/2 
n L+1 , L =p/2+p/(2L) 

TTi+l.i+l = (l-p) 

for the interaction between the pore and its neighboring 
site. 

We repeat, now, the same computation as in the 
sodium case. We first exploit the definition of the sta- 
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tionary measure to write 



:<?2 



1 1 

2?i +2*3 



IL P<lL 
lL ML 



iz pqL+i 



QL-l = \q L -2 + \q L + -^pqL + (b?) 

1 1 1 I. . 

<?l = Tto-i + JP9L + -^Plh + 2 1 1 _ P) r lL 



1 



1 



-^pqL+i + 2~£ ML+1 



= - + (i -p)qL+i 

As for the sodium we obtain 



<ZL-t = <?L 



2iL - i(i + 1) 2iL - i(i + 1) 



41 



21 



-PlL + l 



(B8) 

for £ = 1, — 1, so that the probabilities q\, ... , 

are all expressed in terms of gz, and <zl+i- By combining 
with the expression of qL in (|B7|) we get 



9i = 



2j; 



(l_p)(l-r) 
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FIG. 16. Potassium flux as function of the low-affinity state 
probability in the one-dimensional case for L — 100. The 
points *, x, and + are the Monte Carlo results respectively 
for N K = 3000, 1000, 100. The three lines are the graphs of 
the function (|B2[l with qa given by (|B10|) and r given by (|B9[) 
for the corresponding values of Nk- 



which is obviously equivalent to the last equation in (|B7[) . 

As discussed above, we finally assume qL+i = t/Nk 
and, by exploiting the normalization condition q± + ■ ■ ■ + 
qL + qi+i = 1, we get 



= N K + B + C - y/(N K + B + C) 2 - 4N K C 



and 



qt 



2C 



2p 



(l-p)(l-r) A K 

where we have set 

21? - 3L + 1 2p ( 1 

A = ,B 



(B10) 



1 _ py L + - A A P ),C = l + -pA 



The one-dimensional potassium flux can be finally 
computed by using (|B2|) with qk given by (|B10|) . Monte 
Carlo results are compared with this theoretical predic- 
tion of r) in figure [TBI and the matching is perfect. 

In figure [T7] we have compared Monte Carlo and ana- 
lytical results for L = 1000. Graphs have been plotted on 
the same figures in order to show that even in the one- 
dimensional case the selection effect can be appreciated. 
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FIG. 17. Sodium and potassium flux as function of the low- 
affinity state probability in the one-dimensional case for L = 
1000. The points +, ■, and * are the Monte Carlo estimates 
of the sodium flux respectively for iV N a = 1000, 3000, 5000. 
The points □, x , and o are the Monte Carlo estimates of the 
potassium flux respectively for N K = 1000, 3000, 5000. The 
lines are the graphs of the function (|B1[) and (|B2|I for the 
corresponding values of Nyt a and Nk- 
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